On the Time Average of the Autocorrelation 
Function in Hamiltonian Dynamics 



Pavle Saksida and Tomaz Prosen 

Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000 
Ljubljana, Slovenia 

E-mail: Pavle . SaksidaSfmf .uni-lj . si , Tomaz .ProsenOfmf .uni-lj . si 

Abstract. Rigorous lower bound on the time- average of the autocorrelation function 
of an arbitrary observable is proven in terms of conserved quantities and ergodic 
decompositions of the Hamiltonian dynamics. Improvements with respect to the 
bounds given by Mazur and Suzuki are discussed. 

PACS numbers: 02.3G.-f, 05.45.-a 

AMS classification scheme numbers: 37K05, 37A25, 37A60 



On the Time Average of the Autocorrelation Function in Hamiltonian Dynamics 



2 



1. Introduction 

Understanding the physical implications of ergodic and mixing properties of Hamiltonian 
dynamical systems is one of the main issues of non-equilibrium statistical mechanics. 
It is evident that a complete integrability in the sense of Liouville, or existence of even 
a small (non-complete) set of non-trivial global constants of the motion, is sufficient 
for a Hamiltonian system to be non-ergodic and non-mixing. It has been suggested by 
Mazur [1], and later followed up by Suzuki [2j, that the existence of conservation laws 
can also be connected to divergent transport coefficients expressed - due to a linear 
response theory of Green and Kubo (see, e.g. [31 H]) - in terms of the integrated time 
auto-correlation of the current observable. Namely, Green-Kubo formula expresses the 
conductivity ka with respect to a certain current observable A (say, the energy-current 
for the thermal conductivity) as 



where (3 is the inverse temperature, A^ denotes the time evolution of A. More precise 
meaning of the notation will be explained later. 

Mazur has pointed out an inequality between the time-integrated or time- 
averaged auto-correlation of an arbitrary observable (which may also be interpreted 
thermodynamically as an isothermal susceptibility), and an algebraic expression which 
depends solely on the overlaps between the observable in question and the conservation 
laws expressed in terms of phase space integrals. His theorem is essentially a 
straightforward consequence of an ergodic theorem of Khinchine [5] and demonstrates 
that completely integrable systems should generically behave as ideal (ballistic) 
conductors of heat, electricity, etc. 

Nevertheless, as such a result might seem quite natural from the point of view of 
dynamical systems and ergodic theory, it raised a lot of surprise and attention in the 
solid state community [HIE]- Namely, as expressed in the language of solid state physics, 
one-dimensional completely integrable many-particle systems (with a large number of, 
or in the thermodynamic limit, an infinite number of degrees of freedom) generically 
posses a finite Drude weight, thus a divergent zero-frequency (d.c.) conductivity, and 
therefore behave as ballistic conductors, at all temperatures. Ballistic transport in a 
strongly interacting system at a non-vanishing temperature is certainly a statement 
which raises eyebrows of a solid state physicist. However, complete integrability is not 
even necessary to have such a striking physical implication, it is enough that at least 
one non-trivial conservation law exists which overlaps with the current [8]. 

In this paper we present some results which provide sharp and easily computable 
bounds of the autocorrelation function of an observable in terms of the conserved 
quantities of the system. Mazur's theorem can be considered as a linear (or first order) 
case of our bounds. Under certain conditions, which are quite likely to be fulfilled in 
many concrete cases, our results actually yield precise values of the time average of the 
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autocorrelation function. Mazur's proof relies heavily on a rather involved probabilistic 
and measure-theoretic statement of Khinchine |5]. Our method of proof is completely 
different, rather simpler and it avoids deep results of Khinchine entirely. Nevertheless, 
this method seems to give a good handle on the problem. In particular, it explains how 
the ergodic properties of the system on the one hand, and the form of the observable in 
question on the other, affect the sharpness of the bounds. We expect that saturability 
of our bounds can be an interesting test (probe) of complete integrability. We intend to 
address this issue in another paper. 

In order to have geometric concepts well defined, we shall assume that the number 
of degrees of freedom is finite. Our results are essentially valid for the systems 
with the infinite number of degrees of freedom, which could be argued by taking the 
thermodynamic limit. 

Our setup is the following. Let {M,(jj,H) be a Hamiltonian system, where M is 
the phase space with the symplectic structure u and the Hamiltonian H. Let 

A: M — >R 

be an arbitrary observable. By A^ we denote its time evolution given by 

A\m) = A{jUt))-M-^R, 

where jrn{t):M. — M is the solution of the system {M,uj,H) satisfying the initial 
condition 7m (0) = m. 

The object of our interest in this paper will be the time average of the 
autocorrelation function 

CiA)=lim 1 f{A\A')p dt. 
The symbol (— , —)p denotes the inner product on the L^- space L'i{M) given by 




and 

Z{/3) = [ e-^"^""^ dm 

J M 

is the partition function which we shall assume to exist for strictly positive values of 
(3 (that is, for all finite temperatures if the system belongs to the realm of statistical 
physics). The measure dm is given by the top exterior power uj^ of the symplectic form, 
where r is the number of degrees of freedom. In local canonical coordinates: 

dm = dqi ■ ■ ■ dq^dpi ■ ■ ■ dpr- 

The purpose of this paper is to give simple geometric descriptions of C{A) in terms 
of a set of, say k, functionally independent conserved quantities 

H = H^,H2,...,Hk:M 
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Our main result is the following bound on the time average of the autocorrelation 
function: 

c{A)> (A^^..i^^)/.(H-^)I.n(Ai^^■..i^^)/3, (1) 

l=(!l,...,ij.), 0<|l|<d 
n— (ri]^ , .. .,n^ ) , 0< |n| <d 

where |n| = ^^^1=1 ''^i- The elements of the matrix H are given by inner products 

in a suitably chosen order. Immediately after proving the above bound, we identify the 
conditions under which the bound becomes an equality. We note that the order d can 
be infinite. 

Note that the strict positivity of C{A), implied by the strict positivity of the right- 
hand- side of ([1]), implies ballistic transport, if A is the corresponding current observable 
in the Green-Kubo theory. 

The inequality ([1]) is a non-linear improvement of the inequality 

CiA) = hm i / A')f3 dt > Yl ■ (H-^).,, ■ {A, H,)^, (2) 

T Jo f-^^ 

originally given by P. Mazur in [Ij. Above H denotes the matrix with entries 

Hij = {Hi, Hj)p (3) 

and is its inverse. 

In [1] , working in the framework of statistical mechanics, Mazur treats the evolution 
t A(7m(t)) = v4* as a stochastic process. The main tool he uses is the power spectrum 
/(cj) of the process A^. The power spectrum is given by the Fourier transform of the 
correlation function 

that is by 

1 r°° 

I{uj) = — / (j){t) e"**^ dt. 



27r_ 

The essential ingredient of Mazur 's proof is the well-known result of Khinchine ([5]) 
which states that for every function B: M ^ W we have 



C{B) = hm i / {A', A)^ = /(0+) - /(0-) > 0. 



From the above, Mazur's inequality follows almost immediately. The analytical crux of 
the matter is indeed contained in the above deep result of Khinchine. 

Our treatment relies on direct and simple geometric considerations and does not 
resort to the theory of stochastic processes, let alone to the Khinchine's result. Despite 
its simphcity, our geometric approach enables us to get a better grip on the internality. 
We give simple and meaningful conditions for the inequality to be saturated. We also 
show how the ergodic properties of the system {M,u, H) affect the quantity C{A). 
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Throughout this paper the observables A: M — M will be assumed to be elements 
of Lp{M), where in general LFp{M) is the L^-space generated by the measurable real 
functions /: M — > R which satisfy the integrability condition 

j |j|P e-^^^'^^dm <cx). 

Assuming that the partition function Z{l3) is defined, the measure e~^^^"^^dm is finite. 
Therefore, we have the inclusion Vj^{M) C L'^^{M), and thus A will also be an element 
of L2(M). 

We will divide the discussion into two parts. Firstly, in section 2 we shall study 
the so-called ergodically regular case. This means that the dynamics is ergodic on the 
joint level sets of the conserved quantities. There the bound ([T]) is proven in theorem 
[21 whereas in theorem [1] we give a useful expression of C{A) as the L2-norm of a certain 
projection of the observable A. Theorem [3] identifies the situation in which the bound 
([1]) is saturated. It is reasonable to expect that the result of the theorem [3] will be useful 
in the context of the algebraically integrable systems, where the conserved quantities are 
expressible as polynomials or analytic functions. We list some sufficient conditions for 
the saturation in the remarks following theorem [31 Secondly, in section 3 we study the 
general case, without the assumption of ergodic regularity. The central technical result 
here is lemma [H which is needed to prove the bound iQ, (see corollary [2l) . In theorem 
[Hwe prove that C{A) is equal to the L2-norm of the orbital average of the observable A. 
In proposition [2 we show how ergodic decompositions can be used to further improve 
the bound (P). 



2. The ergodically regular case 

Let, as above, {M,uj, H) be a Hamiltonian system with k < 2r = dim(M) functionally 
independent conserved quantities 

H = HuH2,...,Hk: M 

which are not necessarily in involution. Let 

n = iHi,...,Hk): M^DCR" (4) 

m I — > H{m) = (ai, . . . , afc) 

denote the moment map. We shall first consider the systems whose ergodic behaviour 
is simple in the sense that it is completely determined by the conserved quantities. 
For the sake of brevity we shall denote the level sets of the moment map by 

Ca = n-\a). 

Definition 1 The system {M,uj,H) is called ergodically regular, if it is ergodic on level 
sets Ca for almost every a E D. More precisely, for almost every a E D and for almost 
every m E Ca we have 

lim ^ I B{-im{t)) dt = ] [ B{ma) dm„, 
T^oo T Jo Vol(£„) Jc^ 
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where B G L^(M) is an arbitrary observable, and drria is the generalized microcanonical 
measure on induced by the measure dm which is given by . 

A general example of ergodically regular systems are Liouville integrable systems. For 
the definition of the induced microcanonical measure, see Pj. In the above definition 
it is assumed that every volume Vol(£a) is finite, but a moment of thought shows that 
this is a consequence of the finiteness of the partition function 

The key ingredient in the study of ergodically regular systems is the averaging map 

Let Lj^{D) be the L^-space of measurable functions on D with respect to the measure 
fiD given by 

dfiD = ^ Vol(£,) e-^^''^") da = ^ Vol(£„) 6"^"^ da. (6) 

Then the assignment ([5]) defines the operator 
T: Li(M) L^iD) 

(7) 

B ^ B^. 

Since our measures on M and D are both finite, we have L'^^{M) C Lp{M) and 
L^(D) C L}p{D). Obviously, the averaging operator descends to the map 

r : Lj(M) LliD) 

between the L^-spaces. Let now the operator 

n:L%D)-^Vp{M), z=l,2 

be given by 

{n{G)){m) = G{n{m)). 
The composed operator 

V = TZoT : Ll{M) — ^ Lj(M) (8) 
is then clearly a projector. 

Proposition 1 The operator V given by ^ is a continuous orthogonal operator and 
therefore symmetric on L'^{M). 

Proof. To prove the orthogonality, let F be an element of the kernel of V, and let G 
lie in its image. Then 

[ F{m) G{m) e'^^^"') dm = ! da f F{m^) G{m^) g-^^^"^") dm^ 

Jm Jd JCa 

G^{a) e-^^''('^) I F{mo) dm,) da 



D 



0. 
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Considering suitable convergent sequences in L^i^M) it is easily seen that the subspaces 
Ker('P) and Im('P) are closed, which implies that P is a continuous orthogonal projection 
and therefore a symmetric map of L^(M) into itself. 

□ 

From the above proposition we immediately obtain the basic geometric description 
of the time average C{A) in the ergodically regular case. By || — ||d we shall denote the 
norm on L'^p{D) with respect to the measure dfiD introduced above. 

Theorem 1 Let A be an arbitrary measurable and integrable observable of the 
Hamiltonian system {M,uj,H). Then for the time average of its correlation function 
we have 

C{A) = \\V{A)\\l=\\A^\\l, 
where \\ — \\d denotes the L2-norm on L'^fj{D) given by the measure 
Proof. The fact that our system is ergodically regular gives 

hmi r{A\A'),dt= hm^ / ^(7™(0)) ■ ^(7„.(t)) ■ e'^^^™) dm 



I v^(^-(0))- ^ / ^(^-W) ■^"''''^'"^ dt) dm 



1 



-n(m) 

A{m) V{A){m) e'^^^™) dm. 



Using our notation and the fact that "P is a symmetric projection we get the first equality 
of the theorem, 

{A,V{A)), = {A,V\A)), = {r{A),V{A)), = 

Finally, for every pair of functions A, B: M M. we have 

{V{A),V{B))p = ^ [ P(A)(m) P(S)(m) e-^^(™) dm 

= ^ X"^"^ n^nm^) V{B){m^) e-'^^("^") dm« (9) 

= ^ X^""^") ^""^"^ e-'^^iVoK/:^) da 
= {A'',B^)n. 
which proves our second equality. 

□ 

We shall now turn to the analogues of the right-hand side of Mazur's inequality. 
These expressions will give us estimates for C{A). In concrete cases these expressions 
are likely to be simpler to calculate than the norm of AJ^. More importantly, they will 
be useful in the treatment of the general, ergodically irregular case. 
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Let us consider the set of all monomials corresponding to the multi-indices n = 
(^1, . . . , -n-fc) of degree d or less 

H'^'Hl^' ...Hl\ < |n| = m + ^2 + . . . nfc < d 

composed of the conserved quantities H = Hi, H2, . . . , Hk of our system. We shall order 
these monomials by a combination of the ordering by degree and the lexicographical 
ordering, 

if I 111 I < |n2|, then rii < n2 ; i/ |ni| = |n2|, then ordered lexicographically. (10) 
By Vd we shall denote the subspace of L'^{M) spanned by the above monomials, 

Vd = span{/7fii72"' • • • K'; |n| < d} C LI{M). 

The fact that the conserved quantities are functionally independent implies the linear 
independence of the monomials. The basis {H^^ ...H^^} of Vd is not orthonormal, 
therefore we will need the matrix of the inner product (— , — )/3 on corresponding to 
our basis. The elements of the inner product matrix H are given by 

IU\ / Tjh Tjh Tjlk C7"l ZT^-Z lj1^k\ 

[ti)o{\),o(n) - {Hi H2 ...Hf^ , Hi H2 ...Hj^)(3 

= [ i7{i+"i(m)//^^+"2(m) . . . //^^+"^(m) e-^^^""^ dm, 
Jm 

where o(l), o(n) are integers given by the ordering (fTOl) . The matrix H is non-singular 
due to the linear independence of the monomials. 

Theorem 2 Let A be an observable on an ergodically regular Hamiltonian system 
{M,uj,H) with the additional conserved quantities H2, . . -Hk- Then for every positive 
integer d we have 

u{d) 

C{A) > Yl ^^^2^ • • • (H-^)o(i),o(n) {A, H^^nr . . . Hl^)p, (11) 

o(I),o(n)=0 

where v{d) denotes the number of different monomials of degrees ranging between and 
d in k unknowns. 

The bound (fTT!) could be of practical importance, since it is relatively easily calculable 
in many cases. We notice that the left-hand side of ([1]) is independent of d, therefore 
theorem [2] has the following immediate corollary. 

Corollary 1 Let {M,uj, H) and A:M^M.beas above. Then 

00 

C{A) > Yl H['H'2 ■ ■ ■ H[-)p (H-^)o(i),o(n) {A, H^^H^^ . . . H^^),. (12) 

o(l),o(n)=0 

□ 

Clearly, the bound (|T2l) is sharp. From the practical point of view, it is in general less 
useful. The main problem is the evaluation of the inverse of the infinite matrix H. This 
can be tackled by replacing the monomials by products of suitably scaled orthogonal 
polynomials. 
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Proof of Theorem [2]. Let U = P(L^(M)) C Lij^i^M) be the image of the projector 
v. This is a closed subspace in L'^p{M) and the inherited inner product gives it the 
structure of a Hilbert space. Let 

be the orthogonal projection. We shall prove that the expression on the right-hand 
side of (11 II) is equal to the norm of the vector n{V{A)). Since by Pythagoras' theorem 
any orthogonal projection of a vector is shorter than the vector itself, f|TT]) will follow 
immediately from theorem [H 

Let us introduce a shorter notation 

for the elements of the basis of Vd, and let {h*; i = z/((i) + 1, z/((i) + 2, . . .} be the basis 
of the orthogonal complement Vj- C U composed of the vectors from the basis dual to 
{hn}nen- This means 

{hlh,)p = 0, for ^ = 0...,z/(d), j = u{d) + l,u{d) + 2, . . . . 

Let 

j=0 j=v{d)+l 

be the orthogonal decomposition of V{A) with respect to lA = Vd® V^. Taking the 
inner product of this expression with each of hk for /c = 1, . . . , v{d) gives 

v{d) 

Inverting the matrix Hj ^ = {hi, hk)^ yields the expression of the coefficients q: 

u{d) 

c, = ^(H-i),, {V{A),h,)p. (14) 



3=0 

|2 



The inner product matrix of C L^fj{M) with respect to the basis {hi, . . . , h^(^d)} is H, 
thus the squared norm ||7r(P(yl.)||| of the orthogonal projection of V{A) onto Vd is equal 
to 

/ Co 

||7r(P(A)||J = (co,...,c.(,))-H- ; 

V Cu{d) 

By means of (iMl) we can express this quantity in terms of the products {V{A),hi) p. 
The inner product matrix H is replaced by ■ H ■ H^^, so we get 

h{V{A))r = {{V{A), ho),, {V{A), K^d))p) ■ ■ : 

\{V{A),K^d))p. 
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To complete the proof we only have to show that {A, hi)^ = {V{A), hi)^ for every 
observable A on M. But, the monomials /io(n) = H^^H2^ . . . H^^ are constant on every 
level set £„, and so 



D 



{A,h,)p= / A(m)/i,(m) e-^^(™) dm 

J M 

= /" /i]^ (a) e"^^"" A(m„) dm„) da (15) 
A^{a)h^{a) e-'^"^ da 

From ffTOl) and from the fact that Villi) = h^, we now finally get 
{A^,hl^)n = {V{A),V{h,)), = {V{A),h,), = {V{A),HrHr ...K"),, 
which concludes the proof. 

□ 

We shall now clarify the question, when the inequalities in (fTTi) and in (fT2l) are 
saturated. The answer is given by the theorem [3] bellow, which is more or less an 
immediate corollary of theorem [21 

Theorem 3 (A) Suppose an observable A on M satisfies the two equivalent conditions: 

(i) The function D ^ M. is a polynomial of degree d in the variables (ai, . . . , a^). 
(a) The observable A be expressible in the form 

u(d) 

A{m) = J2 Co(n)(m) H^'im) ■ ■ ■ Hl\m), (16) 

o(n)=0 

where 

c^n)(a) = .r }( r \ / Co(n)("^a) dm„ = const, for every n. (17) 

\IO\.\L.a_) JCa 

Then 

u{d) 

C{A) = Yl ■■■H[')p (H-^)o(i),o(n) {A, H^Hr • • • i^? (18) 

o(l),o(n)=0 

(B) If A^:D ^ W is an analytic function, or alternatively, if A is expressible in the 
form 

oo 

A{m)= J2 Co(n)(m) //r'(m) (m), 

o(n)=0 



where the coefficients Co(n) again satisfy the condition [Tl\) , then we have the equality 

oo 

C{A) = ^^^2^ • • • ^^)/^ (H-^)o(i),o(n) {A, H-^H^ . . . i^D/^- (19) 

o(l),o(n)=0 
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Proof. First we check the fact that the conditions 1. and 2. of the theorem are indeed 
equivalent. Let the observable A be expressible as 

u(d') 

A{m) = doUm)H'l^{m)---Hl\m). 

o(n)=0 

For every m, such that T-[{m) = {ai, . . . , a^), we have Hi{m) = ai on Therefore, 

A^{a)= J2 <n)(«)«r---«r- 

o(n)=0 

This function is a polynomial precisely when all d'^i^^-^{a) are polynomials. In such cases 
the above function can be rewritten in the form 

uid) 
o(n)=0 

where the constants Co(n) are coefficients of the polynomials <i^n-)(a)- Clearly, every do{n) 
is of the form do{n){Tn) = ^ Cj{m)H'[^[m) ■ ■ ■ Hl.''{m) for some choice of the multi- indices 
(ri, . . . Tfc), and for every Cj{m) we have = Cj, which proves the equivalence of 1. and 
2. 

To establish the validity of (A), let 

uid) 

v4^(ai, . . . , ttfc) = ^ Co(n) a^^aa" • • • a^' 

o(n)=0 

be a polynomial. Then the pull-back of A'^ on M is given by 

iy{d) v{d) 

(£(A^))(m) = (P(A))(m) = ^ (c„Hifr/fr---^r)M= E Co(n)(m)/i„(„)(m). 

o(n)=0 o(n)=0 

According to (fT3|) this means that 'P(A) lies in the subspace Vd of the space L'^{M), 
and thus 

vr(P(A))=P(A). 
As seen above, we then have 

C{A) = \\V{A)\\l = \\n{V{Aml 

uid) 
o(l),o(n)=0 

Part (B) of the theorem is an easy consequence of part (A). 

□ 

Bellow we collect some remarks and comments that illustrate the condition f|T71) . 
Remarks 1. The coefficients Co(n)(^) in the expression ^E) are, of course, in general 
not constants. They can be rather arbitrary functions that change in the fibre direction 
of the moment map TC: M ^ D , as well as in the directions transversal to the fibres Ca- 
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2. Let a E D C M} be a regular value of the moment map Ti: M ^ D. Then there exists 
a neighbourhood a E U G D such that the open subset lHr^{U) G M is diffeomorphic 
to U X C and C is diffeomorphic to La- Open sets W C Ti^^iU) = U x C can be 
coordinatized as W = {(ai, . . . , a^, (ii, . . . , [3i)} , where (ai, . . . , a^) are coordinates on 
U G D and {(3i, . . . Pi) some local coordinates on a patch of C. Then the conditions 

Ac,(„)=0, i = l,...,k, o(n) = 0, . . . , (20) 

are sufficient for ( [77[ j and hence for [T^) . Clearly, the above conditions are not necessary 
for ( (77| j. The functions Co(n) o^re allowed to vary in the a-directions, only their averages 
over Ca have to be constant. 

3. Let M be an almost Kdhler manifold. This means that it is equipped with a metric 
g{—, — ) compatible with the symplectic form in the sense that there exists an almost 
complex structure J on M such that 

The most common examples of such manifolds are cotangent bundles M = T*N over 
Riemannian manifolds N. The metric qn on N is extended in the natural way onto the 
metric on the tangent bundle TN , and the symplectic form iJ^ on TN is given as the 
exterior derivative u'^ = d9 of the tautological 1-form 

Om{Vm) = {gN)7r{m){i.DmTT)iVm),rnj , Vm E T„,{TN). 

Here n-.TN N is the natural projection. (Note that m is a tangent vector, 
m E TT^{rn)N .) One can use the metric again to pull the form uP- back to the cotangent 
bundle T*N. Then we have the well defined gradient vector fields VHi associated to 
functions Hi, which are g-orthogonal to the Hamiltonian vector fields Xh- Moreover, 
the gradients VHi are orthogonal to the entire fibres Ca- Therefore, the conditions 

(V//,) ico(n)) = 0, i = l,...,k, o(n) = 0, . . . , u{d) 

are equivalent to the conditions ^2^) . Similarly as the conditions ^2^1 . they are sufficient 
but not necessary for (11). 

4. The condition that is a polynomial is also equivalent to the condition 

^/ |n| = m + na + . . . + rifc > rf, then {A, H^'H^^ . . . Hl'')p = 0. 

This follows immediately from the expansion fT^) and from the equivalence of the 
conditions 1. and 2. of theorem\^ In practice, the problem with the above condition is 
that infinitely many integrals have to be checked. But on the other hand, in some cases, 
the integrals of the functions A ■ H^^ ■ ■ ■ H'^'' might be rather easily computable in some 
concrete cases. 



3. The general case 



The ergodicaUy regular systems treated above are very special. The invariant subspaces 
on which a Hamiltonian system is ergodic can in general be rather wild, and even within 
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a single system they can be of very different types, e.g. full level sets £q,, subsets of 
of lower dimensionality, as well as subsets of having even fractal dimensions. 
One should recall for example typical situations of smooth perturbations of integrable 
systems landing in the context of the Kolmogorov-Arnold-Moser theory. In general, it is 
impossible to parameterize invariant ergodic sets by some manageable (say Hausdorff) 
space over which one could integrate. For these reasons we have to modify our approach 
in order to prove our bounds for a general Hamiltonian system. In particular, for a 
useful description of the ergodic decomposition of our system, we shall revert to an 
inverse limit type construction. This will enable us to prove the results analogous to 
those from the previous section, but valid for the general Hamiltonian systems. The 
analogue of theorem [1] is the following: 

Theorem 4 Let {M,uj,H) be an arbitrary Hamiltonian system with a well defined 
partition function and let A: M be an element of Vp{M). Then for the orbital 

average 

T 



A{m) = hm I- [ A{jmit)) dt 



we have 

C{A) = \\A\\l. 

Before giving the proof, we shall describe the ergodic decomposition of the 
Hamiltonian system {M,uj,H) which will be used in the proofs. In our construction 
we shall use the invariant measure on M given by 

As in the previous section, the crux of the proof will be the replacement of the temporal 
averages by the spatial ones in the context without ergodic regularity. To be able to do 
this, we have to decompose the space M into a collection of invariant sets on which our 
system is ergodic. We will construct such ergodic decomposition by means of successive 
approximations. By definition the Hamiltonian system (M, cj, H) is not ergodic, if there 
exists an invariant measurable set C M such that 

< fi{N) < 1, 

where the inequalities have to be strict. Let the first approximation of our ergodic 
decomposition be a finite partition Afi consisting of invariant sets 

NlNl...,Nl^cM 

with the properties 

o<MiV-)<i, KN!nNj) = o, ut,Nl = M. 

In the next stage we decompose each N} into invariant measurable sets 
{A|., A|.^^, . . . , Aj^} with analogous properties. This yields the partition N2 = 
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{Nl, . . . Nl^} which again satisfies the stipulations analogous to those hsted above. We 
continue the procedure and obtain a sequence {N'n\nm of partitions in which every term 
Un = {N^,..., NJ!J satisfies 

< fi{Nl') < 1, ^(iVf n A^) = 0, uti = M. 

Let now {a„}„gN be a sequence of natural numbers such that for every i G N we have 

1 < Ui < ki. 

Then the set 

Nfa } = lim A^" = U°°_i A^" 

is measurable and invariant. Moreover, the system {M,uj, H) is ergodic on every Ns^an}- 
Remark The sets N^^a^} (^i"^ subsets of level sets C^, therefore they are sets of measure 
zero. These sets can be strange, of noninteger Hausdorff dimensions etc... The collection 
of N^a„} is parameterized by some subset of the set of sequences {a„}„eN of the form 
described above. In order to simplify the parameter set, one could make the above 
construction in a "binary manner" by decomposing each set into only two invariant 
sets Nq^^ and A''""'""'^. The sequences {a„}„,gN would then be the maps 

{an}nm : N — > {0,1}. 

Now we construct the sequence of measurable functions 

An-. M — >m 



by the rule 

if meN^ then = — 1— / A{h) dh = C^..^. (21) 

Vol(A^f) i^n 

Proof of Theorem |4l Our strategy here will be analogous to the one used in the proof 
of theorem [H If A: M — R is an element of L^(M), then by Birkhoff's theorem (see, e. 
g. [To]) the function A is also measurable and A G Lp{M). We clearly have 

lim ^ [ I(7,n(t)) dt = Aim), 

which means that the operator 

A: L|(M) L|(M) 

A I — > A 

is a projector. The image of A are the those functions in L'^{M) which are constant on 
every i^-orbit. This is a closed subspace in L^(M), therefore ^ is a continuous operator. 
Moreover, it is also an orthogonal projector. To see this, let B G ker^, and A = A{A) 
be arbitrary elements in the kernel and in the image of A respectively. We claim that 

{B,A)f}= / B{m) A{m) dfi = 0. (22) 
Jm 
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Let {An}nm be the sequence of functions approximating A as described above. By 
construction we have: For every m E M there exists a unique sequence {a„}„gN such 
that m G Ni^a^^ = hm„^oo N^^. Therefore, 

A{m) = hm An{m), 

since our system is ergodic on every A^{a„}- This imphes 

B{m)A{m) = hm B{m)An{m). 

Functions A and B are elements of L^. Let now Sup be the essential supremum of the 
orbital average A on M. The measure on M is finite, so a function Sup is an element 
of Lp{M). We have the inequality 

B{m)An{m) < B{m) Sup 

which holds for almost every m G M. Therefore, by the Lebesgue dominated 
convergence theorem we have 

B{m) A{m) d/i = lim / B{m) A„(m) d//. (23) 

Since A^ takes the constant value Cj^„ on every A^" for i = 1, . . . , A;„, and since all the 
sets A^" as well as the measure fi are if-invariant, the Liouville theorem gives 

B{m) An{m) d/i = Y 'Ci,„ / B{m) d/i 



M 



i=i 



B{^m{t)) d/i (24) 



for every t G M. Since A{B) = B = 0, we get 



B{m) An{m) d/i 



M 



^C, hm \ {\ B{jUt)) d/i) dt = 0, 



which together with fl23l) proves (!22|) . 

Now, every orthogonal projection is a symmetric operator, therefore 

C{A) = hm i( / A{m) A(7^(t)) d/i) dt 

= {A,A{A))p = (^A,A\A))p 
= kA{A\A{A)), 

- II /IIP 

— 11^11/3) 

which completes the proof. 

□ 

In most cases the orbital average A is impossible to calculate. Therefore, theorem 
m almost never provides a good estimate for C(yl.), apart from the fact that it ensures 
C(A) to be nonnegative. Therefore, as we have done in the ergodically regular case, we 
shall construct more easily calculable estimates which will use the specific information 
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about the Hamiltonian system in question. In the general case this information has two 
sources, the conserved quantities and some ergodic decomposition {N'n\nm of the form 
described above. 

In the proof of theorem|l]we used the sequence {^n}neN of functions approximating 
the observable A. We shall now modify the approximating sequence, so that it will take 
into account the conserved quantities Hi, . . . ,Hk of the system as well as the ergodic 
decomposition. We shall mimic the approach from the previous section, but will replace 
the averaging over the level sets Ca with the averaging over the intersections 

Consider the functions BnX- M — M for n G N and i G {1, . . . A;^} associated to the 
ergodic decomposition {A4}nGN, and given by the rule 

Vol(4r , -) -^^^r . ^^"""^ ] m G iV;^ 

BnM) = { ■ (25) 

; otherwise 

Let now the sequence {Cn}nen of functions 
C„ : M — ^ M 

be given by 

Cnim) = Y,Bn,iM- (26) 

1=1 

The functions C„ are approximations of the orbital average of A. In addition, their 
behaviour is similar to that of the projection V{A) in the ergodically regular case. More 
precisely, the functions C„ are refinements of V{A) which take into account the partition 

Recall that 

{V{A)){m) = -/ A{ma) drua, for me C^. 

Yo\{Cni,m)) Jcn,^^ 

Since Ca = U^Il^N^^, for every m e Ca, have 



{V{A)){m) = — -— / A{ma) dm„ 
Vol(i:„) tn'iiV", 



i=i "^v.i 



Lemma 1 The sequence of functions {Cn}nm has the following properties, 
(i) For every n eN we have 

C{Cn) = \\CJl>\\ViA)\\l = \\A^\\l. 
(ii) The sequence {||C'n||^}neN is non- decreasing. 
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(Hi) The sequence {||C„||^}„gN is convergent, and 
C{A) = lim C{Cn) = lim \\Cn\\l. 

?i— >oo n~*oo ^ 

Proof. Ad 1. The definition of tlie time average of the correlation function, and the 
fact that the functions C„ are constant on the i/- invariant sets N^^^ give 

C{Cn) = I Cn{m)(\im ^ f Cn{lrn{t)) dt) d// 

= / da"^ Bn,i{rna)Qim^— / 5„,i(7„^(t)) dtj dm„ 

= / da V / Bl -{ma) dm„. 
Jd Jn2^^ 



On the other hand, we have 



\\Cn\\l = / Cn{mY dfi 

J M 

p p 

= daj^ Bli{m^) dm„, (28) 
Jd Jn-^^^ 

which proves the first equahty in 1. 

To prove the inequality ||C„||^ > ||P(74)||^ we use ( 1271) and the usual procedure for 
calculating variances. This gives 

- \2 Vol(iV",) 



0< ^(aK)-P(A)(m) 



=1 

h 



Vol(£, 



^ Volf A^" ■) 

^ Vol(£„) 

for every m E Ca- Integrating the above inequality over M with respect to the measure 
d/i gives 



p "-n 



Vol(£J 



) d/i > / P(Af(m)d/.= ||P(A)||J. 



Vol(iV^ 

Observing that the expression ^n{''^i) ^yo\{c^) average of a function over Ca 

and is therefore constant on yields 

= / (£c^K)Vol(iV:,)) da. 
'^-'^ i=i 

Since Cn{m.i) = Bn^i{m) for suitable pairs of m and mj, we see from fl28l) that the above 
expression is indeed equal to the norm ||C„|||, which concludes the proof of the point 1. 
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Ad 2. The proof that for every n G N we have ||Cn+i||^ > ||C„||^, is essentially the same 
as the proof of the inequality ||C„||^ > ||'P(v4)||^ just given. Above we partitioned 
the phase space into the disjoint union vJll^^N^ = M of //-invariant subsets with 
positive measures. To prove 2. we have to partition every A^" further into the union 
U^i;. A^j"'''^'* = A^f . The actual calculations here are then precisely the same as in i., 
modulo slightly different notation. 
Ad 3. First we observe that for every m & M we have 

lim Cn{m) = A{m). 

n^oo 

Indeed, for every m there exists a unique a such that m G and a unique sequence 
{a„}nGN such that m G N^a„} = HJ^^A^^^ . Since our system is ergodic on the limit set 
^{an} ~ l™n^oo A"^, and since the orbit jm{t) is contained in £„, we have 

lim Cnim) = hm I — r / A{ma) dm^ 

n^oo n-.oo Vol(A^^„J J^„^^ 

= lim i / Ai^Ut)) dt 
= A{m). 

Let again Sup be the essential supremum of the orbital average A on M. Then clearly 

Cn{rn) < Sup; a. e. on M, 

since the values C{m) are averages taken over larger sets than those that yield A{m). 
As we already mentioned, Sup is an element of L^(M). Thus, the sequence {C^}„gN 
together with its point-wise limit A"^ satisfies the conditions of the Lebesgue dominated 
convergence theorem. Therefore, we have 

lim ||C„||| = lim / C^(m) d/i = / A (m) d/i = 

n-»oo ^ n^ooj^ J 

This, together with the result of theorem HI concludes the proof. 

□ 

An immediate consequence of lemma [1] is the fact that all the inequalities proved in 
the previous section for the ergodically regular systems hold for all Hamiltonian systems 
without restrictions. In particular, we have 

Corollary 2 Let {M,uj,H) be an arbitrary Hamiltonian system with k conserved 
quantities H = Hi, H2, . . . Hk- Then for every observable A: M ^ W and every G N 
we have: 

u(d) 

C{A) > J2 ^1^2 ■ ■ ■ Ht)p (H-^)o(i),o(n) {A, H-^H-^ . . . H-^)p. (29) 

o(l),o(n)=0 

The above inequality also holds for d = 00. 

Proof. We have just proved that for a general system we have C{A) > V{A). 
Proposition then follows from theorem [2l 
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□ 

We notice that the estimate in the above proposition does not reflect in any way how 
far our system is from being ergodic on the level sets The non-ergodicity is reflected 
in the ergodic decomposition {Nn}nmi cind the information given by {N'n}nm can be 
used to improve the bound (|29l) . In our setup it is quite easy to plug the decomposition 
Nn into fl29|) . Recall the definitions fl25|) and (|26|) . Since for i ^ j the supports iV" and 

of the functions i?„^j and B^^j are disjoint, we have 

{Bn,i,Bn,j)p = Q. (30) 

The function C„ = Y^ili ^n. i is a sum of orthogonal vectors, therefore, 

llC'nll^ = ^ ||-Bn,i||^- (31) 
1=1 

As in the previous section, we can project orthogonally the function G L^(M) on the 
subspace Vd C L'^fji^M) spanned by the monomials iff^i^a^ . . . H^'' of degree d or less. 
We have proved that the vector C„ is longer than V{A), therefore for the projections 
by tt: LI{M) Vd we have 

vr(C„) > 7r(P(A)). (32) 

Taking into account fl30|) and fl3Tl) . and lemma [H we get the following proposition. 

Proposition 2 For every n G N we have 

^ E( E (^n,-^l'^2^ • ■■H[')p ■ (H-^)o(i),o(j) ■ {Bn,,Hi^Hi^ . . . Hi^) ^) . 

i=l o(l),o(j)=0 

The inequality also holds for d = oo. 

□ 

The quality of the above estimate increases with increasing n and d, and for every n >2 
the above estimate is better that the estimate in theorem [2J 
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